Theory of the power spectrum of spin-torque nanocontact vortex oscillators 
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Spin-transfer torques in magnetic nanocontacts can lead to self-sustained magnetization oscillations that in- 
volve large-amplitude gyrotropic vortex motion. This dynamics consists of a steady state orbit around the 
nanocontact, which is made possible because the intrinsic magnetic damping is compensated by spin torques. In 
this article, we present an analytical theory of the power spectrum of these oscillations based on a rigid-vortex 
model. The appearance of vortex oscillations in nanocontacts is not associated with a Hopf bifurcation: there is 
no critical current and the only precondition for steady-state oscillations at finite currents is the existence of a 
vortex in the system, in contrast with conventional spin-torque oscillators that involve large-angle magnetization 
precession. The oscillation frequency is found to depend linearly on the applied current and inversely propor- 
tional to the orbital radius. By solving the associated Langevin problem for the vortex dynamics, the lineshape 
and linewidth for the power spectrum are also obtained. Under typical experimental conditions, a Lorentzian 
lineshape with a current-independent linewidth is predicted. Good quantitative agreement between the theory 
and recent experiments is shown. 

PACS numbers: 75.70.Kw. 75.75.-c, 85.75.-d 



I. INTRODUCTION 

A number of novel phenomena in magnetization dynam- 
ics are made possible by spin-transfer torques. These torques 
arise from the transfer of spin angular momentum between 
a spin-polarized electrical current and local magnetic mo- 
ments, 1 ' 2 and are most apparent in nanoscale structures in 
which current densities are large. An interesting example 
of current-driven dynamics are self-sustained oscillations of 
magnetization, whereby the torques imparted on the local mo- 
ments act in part to compensate the intrinsic dissipation asso- 
ciated with the magnetization dynamics. As a result of this 
compensation, the system can attain a self-sustained oscilla- 
tory state in which the moments precess freely. These self- 
sustained oscillations have been observed in a variety of ge- 
ometries and multilayer compositions, V13 and have drawn in- 
terest for both their fundamental nature and potential applica- 
tions as nanoscale radio-frequency oscillators. 

One class of spin-torque driven oscillations involve mag- 
netic solitons, such as domain walls 14 ' 15 and vortices i 6 ' 13 ' 16 ' 17 
Self-sustained oscillations involving vortices were first 
brought to light in an experiment on spin- valve nanopillars. 6 
In this study, the physical dimensions and aspect ratio of 
the active magnetic "free" layer of the spin-valve were cho- 
sen such that a magnetic vortex ground state is favorable. It 
was shown that steady-state vortex oscillations could be in- 
duced by spin torques with relatively high quality factors, a 
dynamical behavior that differs fundamentally from the tran- 
sient or resonant response studied in other works. 1820 Since 
then, other experiments have shown that similar behavior can 
be observed in magnetic tunnel junctions^ 

A more striking example of self-sustained vortex oscilla- 
tions has been observed in magnetic nanocontacts. 911 In this 
geometry, large current densities are attained by channeling 
electron flow through a magnetoresistive multilayer stack via 
a metallic point contact, which is typically tens of nanometers 
in radius. 4 In contrast to nanopillars where the physical ge- 



ometry confines the vortex motion, the magnetic free layer 
in a nanocontact structure is a continuous film with much 
larger lateral dimensions (typically tens of microns), so trans- 
lational invariance for the vortex is restored within the film 
plane in the absence of structural or magnetic defects. A con- 
fining potential in this medium arises, however, from the Zee- 
man energy due to the Oersted- Ampere fields associated with 
the applied current through the nanocontact. If we imagine 
the nanocontact to be an infinitely long cylindrical conductor, 
then the Oersted- Ampere fields generated by the current flow 
follow a circular pattern in the film plane, which shares the 
same cylindrical symmetry as the magnetic vortex. Because 
of this, the Zeeman energy due to these fields acts as a confin- 
ing potential for the vortex around the nanocontact, with the 
potential minimum being centered on the nanocontact (if this 
is perfectly circular). The resulting dynamics in this geome- 
try, as revealed by micromagnetics simulations, ' ' is a steady 
state motion of the vortex that involves a large orbit around 
the nanocontact, as shown schematically in Fig.Q] By virtue 
of such large orbital motion, the free layer magnetization un- 
derneath the contact undergoes full rotations, which in turn 
results in a large time-varying component of the magnetore- 
sistance. The observation of sub-gigahertz GMR variations in 
experimental nanocontact systems has been attributed to such 
vortex dynamics . 9 ! 1 1 i 21 ~ — 

Nanocontact vortex oscillations represent a dynamical state 
associated with the current flow, because spatial confinement 
for the vortex is absent when no current is applied (aside 
from defects that can pin the vortex). Furthermore, the ex- 
istence itself of a vortex is not guaranteed in the absence of 
currents, because the vortex state bears a higher cost in mag- 
netic energy than the uniform ground state. For these reasons, 
vortex oscillations in nanocontacts differ fundamentally from 
their counterparts in nanopillars, where the physical geome- 
try guarantees both the vortex state and a confining potential 
that allow for damped or resonant dynamics. This distinction 
means that existing theories of spin-torque nano-oscillators, 
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FIG. 1. (Color online) (a) Schematic illustration of current-induced 
vortex oscillations in a magnetic nanocontact system. The vortex, 
assumed to be in the free layer (top) of a spin valve stack, executes 
a large-radius orbit around the nanocontact in a potential generated 
by the Oersted-Ampere field associated with electrical current flow 
through the contact, (b) Vortex orbital motion giving rise to a (c) 
time- varying component of the giant magnetoresistance (GMR) volt- 
age signal. 



which describe a current-driven Hopf bifurcation between a 
damped oscillatory state and a limit cycle, 2732 or other the- 
oretical works describing current-driven vortex dynamics in 
magnetic dots, 3338 cannot be applied directly to the nanocon- 
tact system. 

In this article, we present an analytical theory of the power 
spectrum of current-driven vortex oscillations in magnetic 
nanocontacts. While some elements of this theory have been 
reported in a prior publication, 11 a number of experimental 
observations since have required an extension of that earlier 
work. In particular, the revised theory presented here goes 
toward explaining two key experimental facts. First, it is 
now well established that vortex oscillations in in-plane mag- 
netized systems can be nucleated and sustained without any 
applied magnetic fields. 22 ~ 24 This observation cannot be ex- 
plained by the earlier theory in which a spin-polarization com- 
ponent perpendicular to the film plane, which is expected to 
be negligibly small in the absence of applied perpendicular 
fields, is necessary for sustaining oscillations. 11 While this 
component is known to be important for nanopillars^ the 
reproducible nature of the zero-field result for the nanocon- 
tact system suggests a different spin-transfer mechanism is 
predominant for sustaining oscillations. Second, the quality 
factor of the nanocontact vortex oscillations (i.e., the ratio be- 
tween the oscillation frequency and the spectral linewidth) is 
typically an order of magnitude lower than that for nanopil- 
lar excitations. This is not immediately intuitive, since high 
quality factors have been observed for vortex oscillations in 
nanopillars. 16 Finally, we derive an analytical form for the 
Zeeman energy associated with the Oersted- Ampere field gen- 
erated by the applied currents and show that it varies linearly 
with the radial distance for large amplitude oscillations. This 
linear dependence was used in earlier work, but it was only 
found empirically from numerical calculations. 

This paper is organized as follows. In Section II, we present 
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FIG. 2. (Color online) (a) Geometry of the nanocontact in the film 
plane of the free magnetic layer and system of coordinates used. The 
nanocontact makes a circular cross section with the free layer with 
radius a and centered at the origin. (X, Y) denotes the position of the 
vortex core in the film plane, (b) Profile of vortex core for different 
magnetization tilt angles ©o- 



the geometry, model and equations of motion in the limit of 
the rigid vortex approximation used. There, we derive the 
Thiele equation describing the vortex dynamics. In Section 
III, we consider the large amplitude limit of steady-state vor- 
tex motion in which the orbit is far away outside the nanocon- 
tact. In this limit, we derive some simplified equations of mo- 
tion and analytical solutions for the oscillation radius and fre- 
quency. In Section IV, we examine the stochastic dynamics in 
the large orbit limit and derive analytical forms for the spectral 
linewidth. Some discussion of the results obtained and com- 
parison with experimental data and other work are presented 
in Section V A summary of key results is given in Section VI. 



II. MODEL 

The theory we present describes the magnetization dynam- 
ics of the free magnetic layer in a spin-valve stack, as shown 
in Fig. [2^. We do not account for any coupling that may ap- 
pear between the free and reference magnetic layers, and the 
dynamics of the latter are ignored. The geometry we consider 
is defined in Fig.|2^. In our notation, z represents the axis per- 
pendicular to the film plane, a the nanocontact radius, and d 
the free layer thickness. 

The magnetization dynamics is governed by the Landau- 
Lifshitz equation of motion, with an additional phenomeno- 
logical Gilbert damping term and a spin-torque term r st , 



dm „ dm „ 

— + yo m x H eff = am x — - + T sl 
ot at 



(1) 



where m is a unit vector representing the magnetization ori- 
entation, yo = IM)M is me gyromagnetic constant, //rjH e ff = 
-Vm£ is the effective field, and a is the Gilbert damping con- 
stant. There are two main sources of spin torques that are 
relevant in this geometry, r st = Tcpp + Tap (see Fig. [3]). The 
first involves torques related to electron transport perpendicu- 
lar to the film plane (CPP), which is associated with the CPP 
magnetoresistance of the spin valve and affects only the mag- 
netization dynamics within the nanocontact region, 



Tcpp = -<T\1 Pim • p) m x (p x m). 



(2) 
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FIG. 3. (Color online) Sources of spin- torques in the free ferromag- 
netic layer (Fl), with the arrows indicating current flow, (a) Current 
perpendicular-to-plane (CPP) and (b) current in-plane (CIP) torques 
both contribute to vortex dynamics. The hashed regions indicate the 
regions in which the spin-torques act. NC denotes the nanocontact, 
N the metallic spacer, and F2 the ferromagnetic reference layer. 



The unit vector p represents the magnetization orientation of 
the reference layer and the scalar function V accounts for the 
angular dependence of spin transfer, which depends strongly 
on the layer thicknesses and material parameters of the mul- 
tilayer stack. 2,39 / is the applied current representing electron 
flow and 



ft y 
CTl= Pi- ' 



1 



e Mod na 2 



(3) 



is the CPP spin-transfer efficiency, where Pi is the spin polar- 
ization of the CPP current, and M s is the saturation magnetiza- 
tion of the free layer. Under this convention, / > represents 
the flow of spins (electrons) from the free to the reference 
layer along the negative z direction. The second spin torque 
contribution originates from electron transport within the film 
plane (CIP). Such in-plane currents arise in experimental sys- 
tems because the electrical pads that connect to the bottom 
electrode are usually located at a significant lateral distance 
from the nanocontact. Because electrical currents flowing in 
a ferromagnetic metal are also spin-polarized, spin-torques 
arise in regions in which a spatial gradient in magnetization 
occurs. These CIP torques can be separated into an adiabatic 
and a nonadiabatic component, respectively,^^ 



Tap = -(u ■ V)m +/? m x [(u ■ V)m], 



(4) 



which arise from different spin transport regimes relative to 
the spatial magnetization gradients. In this description, the 
spin torques are more conveniently parametrized by an effec- 
tive spin drift current, 



n y . 
u = P 2 -j(r), 

2eM, 



(5) 



where j is the current density of electron flow in the film 
plane and Pi represents the spin polarization of this CIP cur- 
rent.— In order to keep the ensuing calculations tractable a 
simple model for electron flow is used, which follows from 
the schematic illustration presented in Fig. [3j). The applied 
current is taken to flow uniformly into the cylindrical region 
underneath the nanocontact in the free layer, and the subse- 
quent electron flow within the film plane has a uniform den- 
sity normal to the surface of this cylinder. A similar argument 



is applied to the flow within the nanocontact. Within this ap- 
proximation, the current density in the film plane is purely 
radial from the nanocontact, 



/ r/a r < a 

J(r) = r— — x^ 

LTiaa a/r r > a 



(6) 



In what follows, we consider the dynamics of a single mag- 
netic vortex in the free layer. We will not seek to describe 
the nucleation process here as it is beyond the scope of this 
work. (A recent experiment has shown that vortices and vor- 
tex oscillations can be initiated in a reproducible way in such 
nanocontact systems 43 ). Instead, we will focus on how the 
vortex responds under the influence of the spin torques de- 
scribed above. A crucial assumption we make is that the spa- 
tial profile of the vortex remains constant throughout its mo- 
tion. While the vortex core can deform significantly during 
the course of its motion, 44 and that such changes can be ac- 
count for by spin-wave theories or through the inclusion of 
high-order time derivatives in the equations of motion, 45 such 
effects will not be treated here for the sake of simplicity. 

Let (0, <E>) represent the magnetization orientation in po- 
lar coordinates. In the rigid vortex approximation, the spatial 
profile can be parametrized entirely in terms of the vortex core 
position in the film plane, X = (X, Y), for which the origin is 
taken to be the center of the nanocontact. We consider a vortex 
with an angular variation of magnetization in the film plane of 
the form 



0(jc, y; X,Y) = n tan 



\x-Xl 2 



(7) 



where n - ±1 is the topological charge (n — 1 for a vortex, 
n — -1 for an antivortex). The component of magnetization 
perpendicular to the film plane, which describes the vortex 
core profile, is given by a function of the form = 0(||r - 
R||), where r and R denote the position radial vectors in polar 
coordinates, i.e., ||r - R|| 2 = (x - X) 2 + (y — Y) 2 . The exact 
functional form we choose for the vortex core results only in 
small quantitative differences in the resulting dynamics. For 
the sake of consistency for the remainder of this paper, we 
use a modified version of the Usov ansatz for the vortex core 
profile that accounts for the uniform tilt angle, 0o, of the free 
layer magnetization out of the film planer! 



cos©(r;P) 



Z? 2 + (2cos0 o -l)||r-R|| 



b 2 + \ 



■Rll 



0. 



Rll 


<b 


Rll 


>b 




(8) 



where b is the vortex core radius and p = ±1 is the core polar- 
ization, which describes the orientation of the core magneti- 
zation relative to the normal of the film plane. An illustration 
of the core profile for three tilt angles ©o is given in Fig. |2j). 
The tilt angle is important to describe the magnetization pro- 
file under fields applied perpendicularly to the film plane. A 
detailed comparison between simulated core profiles, this and 
other ansatz, are given in Ref.[37l 



The dynamics of the rigid-vortex is derived following the 
method of collective coordinates. 46 We elevate the core po- 
sition to a dynamic variable, X — > X(f), which allows us to 
express the spatial magnetization profile as = 0[x - X(?)], 
with an analogous expression for O. The ensuing dynamics is 
then generated using the spin Lagrangian 



/ 



dV £(&[x - X(t)], <D[x - X(f)]). 



(9) 



By integrating over the Lagrangian density, £,, with the chosen 
magnetization profile for the vortex, we generate the equation 
of motion for each (generalized) coordinate £, in the usual way, 



d dL 8L 



£,ncj 



(10) 



where F^ nc represents the nonconservative forces such as 
damping and spin-torques. 

The conservative part of the Lagrangian dynamics, i.e., the 
left-hand side of ( fTOb , is given by the spin Berry phase term, 
X.b - (Af,/y)<t>(l - cos 0), and the magnetic energy density 
that plays the role of a potential, 14 - 1/(0, <t>). The Euler- 
Lagrange equations become 



„ dX dU 

Gx — + — 

dt dX 



(11) 



The first term on the left hand side is the gyrotropic term, 
where the gyrovector G is defined by 



M. 



I 



<iVsin0(V<DxV0) 



(12) 



and is equal to G = Gz, where G - 2n(p - cos ®o)ndM s /y, 
The second term is obtained directly from the magnetic energy 
density U — J dV If. If we assume that the lateral dimensions 
of the film are much larger than the vortex core, the magneto- 
static energy, in addition to the usual exchange and crystalline 
anisotropy energies, becomes independent of the vortex core 
position. Translational symmetry in the film plane is there- 
fore only broken by the Oersted- Ampere field generated by 
the electrical current applied through the contact, constituting 
the sole contribution to the potential U. For the sake of sim- 
plicity, we assume that this Oersted- Ampere field is equivalent 
to that produced by an infinitely long current-carrying cylin- 
drical conductor, which allows the Zeeman energy associated 
with this field to be written as 



U z = -fi M, 



H 'f 



dVf(r) sin 0(cos cos <J> + sin <p sin <D), 

(13) 

where Hj = \I\/(2na) is the magnitude of the Oersted- Ampere 
field at the edge of the contact, (r, (p) denotes the spatial vari- 
ables (x,y) in polar coordinates. The function f(f) describes 
the variation of the magnitude of the Oersted field as a func- 
tion of radial distance r in the film plane, 



fir) 



\rja r < a 
\a/r r > a 



(14) 



The nonconservative forces we consider are Gilbert damp- 
ing and spin torques, as discussed above. In a Lagrangian for- 
malism, the force due to Gilbert damping Fq can be included 
by means of a Rayleigh dissipation function W, 



f,G 



dW 



aM s d 



fdV [0 2 + sin 2 0> 2 ] , 



which leads to 



F n = -a D 



<-> dX 



dt ' 



(15) 



(16) 



where a is the Gilbert damping constant and D is the damp- 
ing dyadic, 

B = — fdv(ve® V0 + sin 2 0V<D®VO). (17) 

With the core profile chosen, this term is computed to be 

<— > <— » «— > 

D - D I , where I is the 2x2 identity matrix, D - 

nM s d[2 + sin 2 0o ln(L/b)j /y, and L is lateral size of the free 
layer. 

To describe the CPP spin torques, we first identify the rele- 
vant forces using the Landau-Lifshitz equation (HJ and (f2]), 



F©,cpp - 0, 



(18) 



M c 



Fo.cpp = -o-cpp/ — P(m ■ p) sin z 0. (19) 

7 

By applying the chain rule, we find that the CPP torques can 
be divided into two components, Fcpp = Fcpp,h+Fcpp,j_, which 
describe contributions from the in-plane and perpendicular-to- 
plane components, respectively, of the spin polarization unit 
vector p, 



Fcpp.ii - o"cpp^^ 
1. 



)• 



V0 sin O 



(20) 



-V<t>cos20sin<l> 

2 



- CPP. 



j. = -a-cppIP 



— f dV p ± sin 2 VO, (21) 

r Jnc 



where the volume integration is limited to the region of the 
free layer underneath the nanocontact (NC). Without loss of 
generality, we have assumed that the in-plane spin polariza- 
tion component is along the x axis. To account for the CIP 
torques, we note that the adiabatic and nonadiabatic terms can 
be derived by generalizing the time derivatives to convectional 
derivatives using the spin-drift velocity u, in both the Landau- 
Lifshitz and Lagrangian formulations. The adiabatic compo- 
nent can be obtained from the generalized Berry phase term, 47 



<9<D 

~dt 



(1-COS0) 



dt 



+ u- V 0> (1 -cos©), 



(22) 



which leads to a term that resembles the gyrotropic force, 

M C 
F ad = I dV sin0(V(D x V0) x u, (23) 



while the nonadiabatic term requires a proportionality factor 
/3/a, i.e., we make the substitution 



9 d p „ 

TT -» TT + -U- V 



(24) 



for the time derivatives in the Rayleigh dissipation function, 
which reflects the fact that Galilean invariance is not generally 
present for magnetic dissipative processes in real systems. 48 ' 49 
As expected, the resulting form for the nonadiabatic term is 
similar to Gilbert damping, 



F n _ ad = -— jdV (v© ® V0 + sin 2 0VO ® VO) ■ u. 



(25) 



III. STEADY-STATE OSCILLATIONS 

This section is devoted to the description of the large-orbit 
vortex motion in the steady state. It begins with a brief dis- 
cussion on the initial transient dynamics leading to the steady 
state, where it is shown that the parallel component of the 
CPP torques play a major role in driving the vortex out of 
the nanocontact. Next, a simple form for the Zeeman energy 
due to the Oersted-Ampere field is found. This result, com- 
bined with the CIP torques in the large-amplitude limit, leads 
to equations of motion that are solved to give the steady state 
orbital radius and frequency. 

For the initial transient dynamics, the small-amplitude limit 
is considered in which the vortex is assumed to remain within 
the nanocontact area and close to its center. In this limit, the 
CPP torques can be readily evaluated. For the purposes of 
illustrating the qualitative behavior in this limit, it is assumed 
that f — 1 and no external magnetic field is applied so that 
the tilt magnetization angle can be ignored, i.e. ©o = tt/2. 
The parallel component Fcppj] is nonvanishing only within the 
vortex core where the gradient in the polar angle, V0, and 
cos(0) terms are nonvanishing. As such, it suffices to limit 
the integration in (l2"TT i over only the vortex core region, which 
leads to 



Fcpp.ii = <ril—(nbdp) Un2i+ — — ] 



(26) 



This force is independent of the vortex position and acts to 
drive the vortex out of the nanocontact area. A good estimate 
of the perpendicular component of the CPP torques can be ob- 
tained by neglecting the core contribution using the approxi- 
mation sin 2 » 1, Within this approximation, the integral in 
(|2TT i reduces to a simple integral of the quantity V<J>, which 
yields 

Fcpp,, = tril-f [— ) WTT2 (-71 + Xj) + 0(b 2 ). (27) 

In contrast to the parallel component, the perpendicular CPP 
spin torques acts to drive a gyrotropic motion for the vortex, 
with a magnitude that is inversely proportional to the radial 



vortex distance from the nanocontact center. The combina- 
tion of the two CPP spin torques gives a force acting on the 
vortex that leads to a spiraling motion of the vortex out of the 
nanocontact area. 

For the remainder of this section, we will assume that the 
steady-state regime is attained after such transient dynamics. 
We focus on the large amplitude motion of the vortex around 
the nanocontact, which is relevant for describing the power 
spectrum of the steady-state oscillations as observed experi- 
mentally. We will not concern ourselves with the transient dy- 
namics associated with the nucleation process; this is beyond 
the scope of this paper and will be treated elsewhere. 

The Zeeman energy due to the Oersted-Ampere fields can 
be computed without difficulty if contributions from the vor- 
tex core are neglected; we take = 0o in ( [T3l and retain 
only the magnetization variation in <t>. Let (r, <f>) and (R, ip) 
represent the spatial variables (x, y) and (X, Y) in polar coordi- 
nates, respectively. We proceed by separating the integration 
in Eq. [13] into two parts. The first part involves integrating 
over the contact area, which is found to be 



U z ,m = sin(0 
= sin(0 o ) 



pin r*a 

! o) d<p 
Jo Jo 



dr 



r - Rcosi 



yjr 2 - 2rR cos <p + R 2 



[(2a 2 - R 2 )E(if/) + (a + R) 2 K(i//)], 



9a 



(28) 



where if/ = -4aR/(a-R) 2 , mdK(j/f) = £ dz (l-i/rsin 2 z)- 1/2 

and E(\fi) = J dz (1 - t^sin z) 1 ^ 2 are elliptic integrals. The 
second integral, which runs over the region outside the con- 
tact, is infinite. Nevertheless, we can identify the infinite back- 
ground term from the indefinite integral 



Jf*/ 



dr 



a(r - Rcosip) 
Vr 2 - 2rR cos 4> + R 2 

= 4a\r-R\E 



ArR 



(r-R) 2 



(29) 



by noting that E[-4rR/(r — R) 2 ] = n/2 as r — > oo, from which 
it is deduced that 

t/z,out = HoMsHjd sin(0 o ) [2naR + 4a\a - R\E(if/)] , (30) 

with if/ as defined above. In the large amplitude limit the el- 
liptic integrals can be expanded in a power series in terms of 
a/R, which allows the total Zeeman energy Uz - Uz,m + Uz,out 
to be expressed as 



U z ~ iM)M s Hid sm®Q(4naR), 

= (2n M s d sin ©o) \I\R s k(® ) \I\ R. 



(31) 



This linear dependence on the radial distance of the vortex was 
found empirically through numerical calculations in previous 
work, 11 but we provide a solid basis for this functional form 
here. The linear approximation is very good for R/a > 1, 
which can be seen in Fig. [4] where a comparison with the exact 
solution is presented. 




FIG. 4. (Color online) Comparison between the exact form and linear 
approximation of Zeeman energy, (Jz, as a function of scaled radial 
distance R/a, where a is the point contact radius. The inset represents 
a zoom for R/a < 1 . 



For in-plane magnetized systems, such as the CoFe or 
Permalloy films studied in experiment ^ 9 ' 1 1 i 22 ~ — we assert that 
the CPP spin torques can be neglected in the absence of any 
applied magnetic fie Ids. Under these conditions, both in-plane 
and perpendicular components of the CPP torques should van- 
ish for the following reasons. First, in the absence of any 
perpendicular fields, the p ± component should be vanishingly 
small because the magnetization of the reference layer, which 
is either a hard ferromagnetic material or is exchange biased 
by an antiferromagnet, lies entirely in the film plane. This per- 
pendicular component is important, however, for dynamics in 
large perpendicular applied fields. 1113 Second, if the vortex 
core is sufficiently far from the nanocontact, the magnetiza- 
tion gradient V0 vanishes and cos 20 =* cos n — within the 
nanocontact region. As a result the parallel component also 
vanishes, which can be seen by inspecting Eq.l2"TI 

We are therefore led to the conclusion that CIP spin torques 
are the dominant mechanism for self-sustained vortex oscilla- 
tions in the large orbit limit under zero or low applied fields. 
Let cr 2 represent the spin-torque efficiency for the CIP com- 
ponent in this limit, 



02 = P 2 - 



y 



l 



e M s d Ana 2 ' 



(32) 



such that the spin-current drift velocity can be expressed as 
(r > a), 



u(r) = cr 2 I — r. 
r 



(33) 



With this definition, cr 2 I retains units of angular frequency, in 
accordance with the convention adopted in spin torque oscil- 
lation theory. By substituting this functional dependence for 
the adiabatic torques into d23l above and integrating over the 
region outside the nanocontact, we find the simple and appeal- 
ing result 

F ad = G x u(X). (34) 

under the same as- 



Similarly, the nonadiabatic torques 
sumptions are found to be 



F n _ ad =/?Du(X). 




^ad 



\ V 



n-ad 



D 



'Fq 



FIG. 5. (Color online) Schematic diagram illustrating the different 
force terms in Eq.f36] (a) Static forces: Zeeman potential (Fq), adi- 
abatic (F a( i), and non-adiabatic CF„_ a d) torques, for / > 0. (b) Gilbert 
damping, Fq, associated with the vortex motion v. 



The full equation of motion for the current-driven vortex in 
the large amplitude limit is therefore given by 



Gx 
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(36) 



The role of each term is shown schematically in Fig. The 
conic structure of the Zeeman potential leads to a force that 
pulls the vortex towards the nanocontact center, while the 
nonadiabatic torque counteracts this force pushing the vor- 
tex radially outwards. The adiabatic torque acts as a "boost" 
along the circular trajectory of the vortex orbit, while Gilbert 
damping acts to oppose this motion. 

We are interested in the orbital motion of the vortex around 
the nanocontact, so it is more convenient to work with the 
equations of motion for the radial (R) and angular (ip) vari- 
ables, which can be obtained from d36l l by constructing the ap- 
propriate combinations of the coupled differential equations. 
We find 
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(38) 



The stationary solution is obtained by setting the time varia- 
tion in the radial variable to zero, i.e., d t R = 0. From this 
condition, a solution to the stationary orbit radius R(f) - Rq is 
found, 
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(39) 
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(35) 



There are a number of important points worth noting here. 
First, a physical solution for the radius of the steady-state or- 
bit, Rq > 0, exists only for one current polarity, / > 0, which 
corresponds to the case where electrons flow outward from the 
nanocontact. Intuitively, we can understand this as an outward 
"pressure" exerted by the spin torques that counterbalances 
the Gilbert damping and therefore prevents the vortex spiral- 
ing into the contact center, which corresponds to the position 
of energy minimum in the Zeeman energy. Second, the orbital 
radius is independent of the applied current and the nanocon- 
tact size, and depends only on material parameters. As such, it 



also follows that there is no critical current for self-sustained 
oscillations, in stark contrast to nanopillar spin-torque oscil- 
lators that involve large amplitude spin waves. The present 
theory therefore predicts the existence of vortex oscillations 
for any value of the applied electron current / > 0, provided a 
vortex is already present in the nanocontact system. The on- 
set of oscillations in this scenario is therefore be determined 
by a threshold for vortex nucleation, rather than a supercriti- 
cal Hopf bifurcation as in the case for nanopillar oscillators. 
Indeed, the experimental observation that vortex oscillations 
persist below an onset current, after having been nucleated 
above this current, 9,25 supports this idea. 

From the stationary condition for the radial dynamics, the 
frequency of the vortex oscillation can be evaluated directly 
by substituting Rq for/? in (l38T l. If we assume that the nonadi- 
abatic torque is similar in magnitude to the adiabatic torque, 
f3 =a a, which is found to be a good approximation for transi- 
tional metal ferromagnets, 50,51 the first term on the right hand 
side of d38l l can be neglected and the expression for the fre- 
quency reduces to the simpler form 
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FIG. 6. (Color online) Variation of the GMR voltage signal as a 
function of vortex radial distance, in normalized units. The variation, 
computed by numerical integration, is compared with the analytical 
result and simple asymptotic forms for small (|v|/Vo - R/a) and large 
amplitude (Eq,l45ll orbits. 



generality, it suffices to consider the x component of magneti- 
zation in the free layer, 
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where \I\/Ro plays the role of an effective magnetic field, with 
the sense of rotation is given by the vortex core polarization p. 
This is also an appealing result because it allows a quantita- 
tive estimate of the vortex orbital radius to be obtained readily 
from experiment: By measuring the slope of the frequency 
versus current curve in zero applied field (@o = tt/2), a mea- 
sure of Ro can be obtained directly. 



IV. POWER SPECTRUM 

The power spectrum measured in experiment is related 
to voltage variations v(f) associated with the time-varying 
change in magnetoresistance in the nanocontact region. There 
are two contributions to the magnetoresistance: one from CPP 
currents, the other from CIP currents. Because the vortex 
orbital motion leads to the largest time-varying variation in 
magnetization in the contact region, the larger contribution 
to the total magnetoresistance is expected to come from the 
CPP component. Let AV represent the total CPP magnetore- 
sistance for the spin valve. The time-varying component of 
interest can be written as 



yV 2 - 2rR cos((p - ip) + R 2 

(42) 
This integral can be solved using a power-series expansion in 
the variable r/R < 1 and summing over all terms. We find 



v(f) = --na d (AV) sinyp 2 F\ 
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where 2 F\ is a hypergeometric function which can expressed 
by the integral 



2Fi(a;b;c;z) = 
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A plot of the variation of the magnetoresistance signal magni- 
tude as a function of orbit radius is shown in Fig. [6] For small 
amplitude motion, the GMR signal exhibits a linear variation 
as a function in the radial distance. In the limit of large am- 
plitude motion, the hypergeometric function can be approxi- 
mated by 2 F\ \-\\ j\2;(a/R) 2 ] =* 1 - a 2 /(8fl 2 ), which leads 
to a simple expression of the time-varying magnetoresistance 
signal 
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It is assumed that the fixed layer magnetization remains static 
and is uncoupled to the vortex dynamics in the free layer to 
simplify the calculations, although in practice this may need 
to be accounted for. To simplify the calculations, it will also 
be assumed that the fixed layer magnetization is uniform in the 
nanocontact region. As such, the voltage calculation amounts 
to averaging over the in-plane component of magnetization of 
the free layer over the point-contact region. Without loss of 



where Vo = na 2 dAV/2. This result indicates that the vortex 
nanocontact oscillator can be considered as a phase oscillator 
to a very good approximation, with radial fluctuations that de- 
crease like l/R 2 . As such the power spectrum of the voltage 
fluctuations S (w), 



S(co) 



-f 



dt K(f)e 



(46) 



which is the Fourier transform of the voltage autocorrelation 
function K(t) = (v(f)v*(0)), is well described by the phase 



variance of the oscillator, Atp 1 , 

K(t) = V 2 exp [i(<p(t)}] exp 



--A<p l (t) 



(47) 



where A(p 2 (t) = ((p(t) 2 ) - {ip(t)) 2 . Once the phase variance is 
known, the lineshape and linewidth of the power spectrum can 
be obtained directly. 

The power spectrum of any oscillator is broadened at fi- 
nite temperatures by thermal noise. We can describe the influ- 
ence of such stochastic processes on the oscillator dynamics 
by including additional noise terms in the equations of motion 
(|36T l, which we can express symbolically as the set of coupled 
Langevin equations 52,53 



dX 
dt 



V(X) + y/q 77(f). 



(48) 



Note that the noise contribution enters as an additive term, 
rather than a multiplicative process (such as a random field in 
the Landau-Lifshitz equation); it has been shown elsewhere 
that such an additive noise is adequate for of vortex dynam- 
j CSj 52 i 53 ^ _ (jj x ^ jj^ j s a two-component vector that represents 
a Gaussian white-noise forcing, which possesses the spectral 
properties 



(77,(0) = 0; <i7,(f)i7/0> = 25ijS(t - t'\ 



(49) 



with q representing the thermal noise amplitude. The choice 
of 
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(50) 



ensures that the fluctuation-dissipation theorem is satisfied in 
the absence of spin-transfer torques. 

Since we are interested in describing the fluctuations about 
the steady-state orbital motion, we linearize the Langevin 
equations about the stationary orbit Ro, R(t) - Rq + r(t), which 
allows us to obtain from d48l . 



dr 
dt 



~~rT 



r + q\\ 



r 



+ yfqmit), (51) 



where 77^ possess the same spectral properties as 77x,y. By 
analogy with spin torque oscillator theory, we identify a 
restoration rate T r - <T2la 2 /R^, which describes the rate at 
which fluctuations in the orbital radius are damped out. The 
radial fluctuations are independent of the phase fluctuations 
and are driven by white thermal noise, subject to a spurious 
drift. By using material parameters relevant for typical ex- 
periments, it is straightforward to show that the term in q is 
small compared with the additive noise term proportional to 
yfq, which allows the Langevin equation in r(t) to be reduced 
to the simpler Ornstein-Uhlenbeck process, 



dr 

— = -T>+ yfqrf R {t). 

at 



(52) 



The formal solution to this equation, r(f) = cexp(-F r t) + 
y/q exp(-r r ?) J t/T77 S (r)exp(r r r), where c is a constant, al- 
lows the two-time autocorrelation function for the radial fluc- 
tuation to be computed directly. In the limit of large times 



(at which the initial correlations are forgotten), we recover the 
known result 



<r(fM?')> = f exp(TY|f-r'|). 

1 r 



(53) 



The linearized Langevin equation for the phase dynamics is 



dip k\I\ I r(f)\ ,_ 

m = GR- l l -^ + ^ T >* it) ' 



(54) 



where we have neglected a cross term in r(t)i]^(t). In contrast 
to the radial dynamics, which is uncoupled from the phase 
fluctuations, the phase variable is coupled to the radial fluctu- 
ations. As a result, the phase dynamics is driven by the radial 
fluctuations r(t), which appears as colored noise, in additional 
to the additive white noise proportional y/qrj^it). From the 
formal solution to this differential equation, and by using the 
spectral noise properties for 77, we obtain the phase variance 
to be 
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We can identify a nonlinearity parameter v, 
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(55) 



(56) 



along with a "linear" linewidth parameter Acoq = q/R^, which 
allows the expression for the phase variance to be simplified 
to the form, 



1 ? 

-A(p(ty = Awo 
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r, 



(57) 



This expression is identical to the phase variance for a spin 
torque oscillator, which describes an inhomogeneous broad- 
ening of the spectral line due to radial (amplitude) fluctua- 
tions. 32 

Two limiting cases appear as a result of this inhomoge- 
neous broadening. In the "low" temperature limit in which 
the coherence time, r c , of oscillations is much longer than 
the inverse of the restoration rate, T,t c » 1, the exponential 
term can be neglected and the phase variance is proportional 
to \t\. As such, the power spectrum in this limit is described 
by a Lorentzian lineshape with a full width at half maximum 
(FWHM) of 



Awr 



= Aojq (l 



- 2 )=^ 



v 2 )- 



(58) 



Through the linear dependence of Awo on the noise parameter 
q, the linewidth in this limit varies linearly as a function of 
temperature and is inversely proportional to the square of the 
steady-state radius R^. Because the orbital radius is indepen- 
dent of the applied current, the linewidth is also independent 
of the current in the low temperature limit, which is in stark 
contrast to conventional spin torque oscillators for which the 
supercriticality, which describes the ratio between the applied 
and threshold currents, is a primordial factor. For the present 
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case the linewidth is determined primarily by material param- 
eters. In the opposite "high" temperature limit in which the 
coherence time is much shorter than the inverse of the restora- 
tion rate, T r r c <Z 1, the exponential function in the phase vari- 
ance can be expanded in a power series to give, 



^(0 2 -Aw l|f| + ir>-; 
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(59) 



If the nonlinearity v is sufficiently large, the linear term in \t\ 
can be neglected and the power spectrum is described by a 
Gaussian lineshape with a FWHM of 



Aw HT = 2 V21n2 |v| y/Ato^ r ~ 2.35 |v|a ^ cr2/ 



(60) 



In this limit, the temperature dependence of the linewidth is of 
the form T 1 ^ 2 , which is consistent with inhomogeneous broad- 
ening due to phase-amplitude coupling. 32 Furthermore, the 
linewidth also acquires a square-root dependence on the ap- 
plied current through the restoration rate, which is markedly 
different from the low temperature case. This gives a direct 
experimental means of identifying the temperature regime by 
measuring the current dependence of the linewidth. 



V. DISCUSSION AND COMPARISON TO 
EXPERIMENTAL DATA 

The fundamental premise of this theory is that sub-GHz 
voltage oscillations observed in experimental nanocontact 
systems are due to the orbital motion of a single magnetic 
vortex. The quasi-linear current dependence of the oscillation 
frequency provides the strongest evidence to support this hy- 
pothesis to date. However, there remain open questions con- 
cerning how the single vortex state is attained from an initially 
uniform magnetic ground state. This crucial issue is a prob- 
lem of topology and can be discussed in terms of the Skyrmion 
number Q? A A vortex or an antivortex possesses a half-integer 
Skyrmion charge of Q — np/2, where n is the winding num- 
ber, with n — 1 for vortices and n--\ for antivortices, and p 
is the core polarization. As such, a system with a single vor- 
tex possesses a charge of Q - 1 /2, while the uniform state has 
Q - 0. 55 In this light, the single vortex and uniform magnetic 
states are in different topological sectors, which means some 
physical process that does not conserve topological charge 
needs to occur during the nucleation. It is interesting to note 
that a vortex-antivortex pair with parallel core polarizations 
is in the same topological sector as the uniform state, while 
a vortex-antivortex pair with opposite core polarizations pos- 
sesses Q = ±1. The latter case is interesting because it ad- 
mits a rotating solution^^S in the absence of damping, such 
a vortex-antivortex pair (or "dipole") rotates about its center- 
of-mass with a frequency of 0J4 - 4yA/(MJ 2 ), where / is the 
separation between the vortices and A is the exchange stiff- 
ness. Indeed, the possible existence of such rotating vortex 
dipoles in a nanocontact geometry were seen in a recent nu- 
merical study by Berkov and Gorn. 58 While such a solution 
is theoretically appealing because a uniform magnetic state 



is conserved far from the nanocontact, it is not immediately 
apparent how dynamics involving vortex dynamics can ac- 
count for experimental observations. First, there is no obvi- 
ous mechanism by which a linear current dependence of the 
oscillation frequency can be obtained. As we have shown in 
a recent study, 43 the Oersted- Ampere field gives rise to a Zee- 
man potential that drives the vortex pair apart, which would 
lead to a decrease in the oscillation frequency with current. 
Second, a rotating dipole pair possesses a Skyrmion charge of 
Q = ±1, which is also topologically distinct from the uniform 
ground state. Therefore, as for the single vortex case, some 
charge non-conserving process needs to take place. 

As discussed briefly in Section [Hi] we assert that in-plane 
spin-torques, rather than perpendicular-to-plane torques, 
should play a dominant role in the vortex dynamics in mag- 
netic nanocontact systems under zero or low-applied fields. 
This hypothesis is strongly supported by a number of re- 
cent experiments in which vortex oscillations have been re- 
producibly initiated in the absence of any applied magnetic 
field i 22 ! 23 ' 43 Under such conditions the magnetization orienta- 
tion of the reference magnetic layer should lie entirely in the 
film plane, so there is no reason to expect a large contribution 
from the Fqpp,± term, which is necessary for the existence of 
self-sustained oscillations. 11 This point is well-illustrated by a 
recent experiment on vortex oscillations in magnetic nanopil- 
lars in which large out-of -plane fields are required. 13 While 
the dominant CIP adiabatic torques lead to the same functional 
form as the Fcpp,± term in the reduced equations of motion, 
there is a crucial difference between the two: CIP torques act 
independently of the vortex core polarization, while the ex- 
istence of vortex oscillations under CPP torques depend on 
relative orientation of the core polarization with respect to the 
perpendicular spin torque component. 1113 Therefore, vortex 
core reversal would not restrict self-sustained oscillations un- 
der CIP torques but would lead to a damped oscillatory regime 
under purely CPP torques. In this light, it would be interest- 
ing to see whether such core reversal processes could be de- 
tectable in a time-resolved experiment. 

The hypothesis of self-sustained oscillations driven by CIP 
torques might explain the ubiquitous presence of higher har- 
monics in the power spectrum. In this picture, the magnitude 
of the spin-torques due to lateral currents in the film plane de- 
termines the shape of the vortex orbit. A uniform radial flow 
has been assumed in the present work as a matter of simplicity, 
but the actual current distribution in an experimental system 
would certainly be more complex. As such, one could expect 
elliptical orbits to result, which would lead to significant con- 
tributions to the harmonic content of the power spectrum. 

An important test of the model described here is the slope 
of the frequency versus current relation, 



dcj k K 2 aD 
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(61) 



which, as discussed previously, is a material-dependent pa- 
rameter. It is a quantity that can readily be extracted from 
experimental measurements and represents a robust charac- 
terization of the oscillator properties because the oscillation 
frequency, as opposed to other spectral parameters such as 
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the linewidth or power, is a stable physical parameter that 
depends solely on the vortex dynamics and is independent 
of thermal fluctuations. By using the experimentally deter- 
mined values of /vqMj = 1.56 T, d — 3.5 nm, a - 0.013, 
L — 10 /im from Ref. |23, and by assuming a spin polariza- 
tion of P — 0.5 and b — 10 nm, we find a theoretical value 
of doj/d\I\ = 4.4 MHz/mA, which is within a factor of two 
of the observed slope of 7.4 MHz/mAi22i22 For the low-field 
vortex oscillation studies of Refs. |9| and |2U a frequency slope 
of as 30 MHz/mA is predicted (by assuming hqM s - 0.8 T, 
d = 5 nm, a = 0.01, L = 10, P = 0.5, and fo = 10 nm). 
This is a factor of three larger than the observed slope of « 
10 MHz/mA. While some uncertainty exists in the spin polar- 
ization P of the applied currents and the vortex core radius b 
in these experiments, the agreement between theory and ex- 
periment nevertheless remains relative good on a quantitative 
level for a simple model with no adjustable parameters. 

Another important point of comparison is the spectral 
linewidth. Indeed, understanding the underlying physics gov- 
erning linewidths is crucial for any potential oscillator appli- 
cations. From the experiment described in Ref. |22|, it was 
shown that the autocorrelation function for the measured high- 
frequency voltage signals can be well-described by a decaying 
exponential function, (v(f)v(O)) oc exp(-|f|/r c ), where t c is 
the characteristic coherence time of the voltage oscillations. 
As discussed in the previous section, this form for the au- 
tocorrelation function indicates that the spectral lineshape is 
Lorentzian, with a FWHM given by Aw = 2/t c . For an ap- 
plied current of 18.7 mA in zero applied field (@o = n/2), 
a coherence time of 140 ns was observed. 22 From the results 
obtained for the low-temperature limit in the previous section, 
the present theory predicts the current-independent value of 
t c » 70 ns with the same experimental parameters, which is in 
reasonable agreement with the experimental value. In a simi- 
lar manner, a current-independent linewidth of A/ » 2 MHz 
is predicted for the experiment of Ref. 21, which is within a 
factor of three of the experimentally measured FWHM of 0.78 
MHz for the oscillation mode at 128 MHz. The time-domain 
analysis performed for this particular experiment shows that 
the dominant contribution to the spectral linewidth is due to 
phase noise, so the present theory is applicable. 

The discrepancies between the theoretical and experimental 
values of the frequency and linewidth might originate from the 
way the Oersted- Ampere field and spin torques are computed 
in this theory. For instance, the spatial profile of the Oersted- 
Ampere field has been computed by approximating the cur- 
rent flow through the nanocontact with the flow through an 
infinite cylindrical wire. On the other hand, it has also been 
assumed that a large component of this current flows later- 
ally from the nanocontact because of how the contact pads 
are located. These assumptions are not compatible with one 
another, and as a consequence both the spin-torques and the 
Oersted- Ampere field are overestimated in the theory. A sim- 
ple way of correcting these estimates is to include ad-hoc cor- 
rection factors Ko < 1 and cr^ < 1 into the definitions of the 
Zeeman energy and spin-torque parameters, respectively, 



Following this line of reasoning, the frequency versus current 
slope and Lorentzian linewidth also acquire these correction 
factors, 
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but with different functional forms. By setting the ratios be- 
tween the theoretical and experimental values of the frequency 
and linewidth to unity, it is found empirically that kq « 0.4 
and (To w 0.3. The empirical value of kq suggests that a half- 
infinite cylindrical wire, for which one would expect a correc- 
tion factor of 0.5, gives a better approximation for comput- 
ing the Oersted-Ampere field. It would be interesting to see 
whether more detailed finite-element calculations of the cur- 
rent flow in realistic geometries would give such correction 
factors. 



VI. SUMMARY 

In summary, a theory of the power spectrum of current- 
driven vortex oscillations in magnetic nanocontacts has been 
presented. The theory is based on a rigid-vortex model and 
the equations of motion describing the vortex dynamics have 
been derived and solved in the steady-state limit. In contrast to 
conventional spin torque nano-oscillators that involve large- 
angle magnetization precession, the self-oscillatory state in 
the nanocontact system is found to exist with the only con- 
dition being the existence of a vortex in the system. As such, 
the onset of oscillations does not involve a Hopf bifurcation 
and therefore no critical current is predicted. It is found that 
spin-torques due to current flow in the plane of the free mag- 
netic layer are crucial for the existence of self-sustained oscil- 
lations. The oscillation frequency is found to vary linearly as 
a function of current, in accordance with experimental obser- 
vations, with the function form 
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where the orbital radius Rq is current-independent and de- 
pends only on material parameters, 

G 2 cr2a 2 

Ro = ^-sgn(7). 

KaD 

In the low-temperature limit in which the coherence time of 
the oscillations is greater than the restoration rate of the radial 
fluctuations, it is found that the power-spectrum is described 
by a Lorentzian lineshape with a linewidth (FWHM) of 



Awlt = Awo (1 + v ), 



where Ao»o is a "linear" linewidth that is proportional to the 
temperature and visa "nonlinearity" parameter that represents 
the ratio between the magnitude of the damping dyadic and 
the gyrovector, v = aD/G. At higher temperatures, the line- 
shape is inhomogeneously broadened by fluctuations in the 
orbital radius and a Gaussian lineshape is predicted, with a 
current-dependent linewidth that varies like V7. 
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